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Abstract. We discuss bispectra of single-field inflationary models described by general 
Lorentz invariant Lagrangians that are at most first order in field derivatives, including the 
fast-roll models investigated by Noller and Magueijo. Based on a factor analysis, we identify 
the least correlated basic contributions to the general shape and show quantitatively which 
templates provide a good approximation. We compute how relative contributions of basic 
shapes to the total bispectrum scale as slow roll is relaxed. To enable future comparison with 
CMB observations, we provide a modal expansion of these non-separable bispectra in Fourier 
space, employing the formalism by Shellard et al. Convergence is rapid, usually better than 
ninety-five percent with less than thirty modes, due to the smoothness of these primordial 
shapes. 

Truncated polynomial modal expansions have restrictions, which we highlight using an 
example with slow convergence. The particular shape originates from particle production 
during inflation (common in trapped inflation) and entails both localized and oscillatory 
features. We show that this shape can be recovered efficiently using a Fourier basis and 
outline the prospect of future model parameter extraction and N-body simulations based on 
modal techniques. 
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1 Introduction 

Observations of primordial curvature perturbations by measurements of the cosmic microwave 
background radiation (CMB) and large-scale structure (LSS) are in good agreement with 
predictions of inflationary scenarios. However, the number of different models that predict 
almost scale-invariant, Gaussian fluctuations is vast. Non-Gaussianities (NG) allow for fur- 
ther discrimination if they are observed in current or upcoming CMB [1, 2] and LSS [3, 4] 
measurements (i.e. via the scale dependent bias [5-8]). 

Up until recently, the main focus has been on predicting and measuring the amplitude 
of three separable shape templates of the bispectrum (see [9] for a collection of reviews): the 
local one [10-12], easily recovered by means of the (non-linear) 5A^-formalism [13, 14], the 
equilateral one, often used to approximate the bispectrum of DBI inflation [15] or scenarios 
with more general non-canonical kinetic terms [16], and an orthogonal shape which is hard 
to generate directly (see however [17] for a possible source). 

Any detection of primordial NG would rule out simple single-field models, since the 
bispectrum is slow-roll suppressed [18]. The magnitude of the three templates in the data 
can be expressed by a single, scalar amplitude, the non-linearity parameter /nl- A recent 
CMB analysis based on the WMAP7 data set alone [19] yields at the 95% confidence level 

-10 < /^T' < 74, -254 < /^f < 306, -410 < /°f °§ < 6, (1.1) 

consistent with a Gaussian spectrum. The data analysis, using estimators [20] that are 
optimal even in the presence of foreground contamination, depends crucially on separability 
of the bispectrum to reduce the dimensionality of integrals. Such an analysis is not feasible 
for non-separable shapes, since computing time of an optimal estimator scales with O (/max) 5 
here / is the multipole number (/^ax ^ 2000 for PLANCK [1,2]). 

Albeit convenient, relying on templates is problematic. First, a quantitative under- 
standing of how well the three templates approximate real bispectra is needed. In this paper 
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we carry out this analysis for general single-field models (Lorentz invariant Lagrangians at 
most first order in field derivatives). These models were first discussed in [21] under the 
assumptions of slow roll and small deviations from unity of the speed of sound. In [16, 22] 
an arbitrary speed of sound was allowed and in [23] slow roll was relaxed to yield the general 
bispectrum investigated further in this paper. ^ After dissecting the resulting general shape 
into a handful of different contributions/subclasses, we perform a factor analysis to identify 
the least correlated basic shape constituents; we compute the overlap with templates and 
show that a restriction to common ones is usually sufficient for these models (see [24] for 
related work on single-field Galileon models that appeared while completing this article; our 
conclusions are in line with [24]). 

Nevertheless, examples of non-separable shapes that show little overlap with the tem- 
plates are known, e.g. [25, 26]. Furthermore, by relying on a single number such as /atl 
to characterize non-Gaussianities, we waste the opportunity to tell individual sources apart 
(primordial/foregrounds/artefacts) since their respective shape function can act as a unique 
fingerprint. 

To ameliorate these shortcomings, Shellard et al. proposed an expansion of the bispec- 
trum into a complete set of separable modes (polynomials) [27]. The expansion coefficients 
of the late-time bispectrum in the CMB data can be calculated using a generalized version 
of the WMAP algorithm for separable shapes, which scales as O (/max)- Several concrete 
models proposed in the literature have been analyzed [27-30] and tested for convergence, 
with positive results. Among these are, in addition to the ones mentioned above, warm 
inflation [31], oscillatory features in the potential leading to resonances [32, 33] and ghost 
inflation [34]. A related mode expansion for the trispectrum was proposed in [35, 36] and an 
expansion into Fourier modes has been proposed in [37] with the aim to enhance convergence 
for oscillatory shapes. 

The formalism appears promising to perform a blind search for late-time NG signals, 
e.g. caused by cosmic strings [27], and can be used to tell apart foregrounds/artefacts from 
primordial signals in the PLANCK data analysis. Current developments entail the estimation 
of NG signals in the large-scale structure [38] and the generation of non- Gaussian initial 
conditions for N-body simulations [39, 40]. 

In this paper we apply the (polynomial) modal expansion to the general single-field 
models described above. We find that naively chosen basic shapes fall into three categories 
that cannot be separated further observationally at present, due to their strong correlation 
with templates. Using a factor analysis, we identify least correlated basic shape constituents 
that offer better observational separability. Focusing on concrete models, we show that the 
three dominant contributions correspond to the three common templates; the remaining 
distinct subclasses are subdominant and hard to detect due to a reduced but still significant 
degree of correlation. In addition, we identify systematics that complicate reconstruction of 
subdominant constituents. 

We further show how the relative contributions of the subclasses to the total bispectrum 
scale as slow roll is relaxed (the relaxation of slow roll requires fine-tuning to retain a nearly 
scale-invariant power spectrum). As a consequence, one could in principle pinpoint the 
concrete model if non-Gaussianities were observed. Particularly, if expansion coefficients 
were extracted from observations, one could perform a Markov chain Monte Carlo (MCMC) 
analysis to scan the parameter space and identify the best fit. To this end, the identification 

^Deviations from slow roll can still be consistent with scale invariance of the power spectrum due to possible 
cancellations; in [23] bispectra of maximally slow-roll violating models (up to e ^ 0.3) were constructed. 
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of well chosen basic shape constituents is crucial to perform this analysis efficiently: if basic 
shapes are chosen properly, the expansion coefficients of a general shape can be obtained 
from a weighted sum of basic shape coefficients (only weights depend on model parameters), 
which need to be computed only once. 

As an aside, we comment on work in progress on applying modal techniques to N-body 
simulations. 

However, truncated modal expansions have restrictions: we investigate a shape with 
both localized and oscillatory features in the bispectrum [25, 26] that has slow convergence 
for a polynomial expansion. This shape originates from backscattering of particles produced 
during inflation from the inflaton condensate [41, 42], a common phenomenon in trapped 
inflation [43-47] or monodromy inflation [48]. However, we find faster convergence for this 
shape using the Fourier expansion of Meerburg [37]. Furthermore, Fourier expansion coef- 
ficients are better suited to identify key frequencies in the bispectrum, leading to improved 
model identification properties. 

The paper is organized as follows: we start with a pedagogical introduction to modal 
decompositions in section 2. We focus on the expansion of the bispectrum in separable 
polynomial modes as introduced by Shellard et al., and explain the considerably involved 
notation in detail. Section 3 contains a review of the generalized single-field models by 
Seery and Lidsey [21] (section 3.1), Chen et al. [16, 22] (section 3.2), as well as Noller 
and Magueijo [23] (section 3.3). Readers familiar with the modal decomposition may want 
to skip section 2, while readers familiar with general single- field models may want to skip 
Section 3. We separate general bispectra of singe-field models into subclasses, identify the 
least correlated basic shape constituents and perform a modal decomposition in section 4. 
Our results are summarized in section 4.4 and related to recent work in section 4.5. We follow 
with an outlook to future work, especially the application of modal techniques to N-body 
simulations, in section 4.6. Fourier modal expansions are discussed in section 5 and applied 
to non-Gaussianities from particle production during inflation. Some lengthy expressions are 
provided in the appendix. 

1.1 Non-Gaussianities, the bispectrum and the shape function 

Non-Gaussianities entail the study of higher order correlation functions beyond the two-point 
function, or power spectrum. In this paper, we focus on the three-point correlation function 
of the comoving curvature perturbation TZ using the modal expansion of [27, 29] . The Fourier 
transform of the thre-point function, the bispectrum of primordial fluctuations B^{ki, k2, ^3), 
is the object of interest; it is defined as 

(7^(fcl)7^(fc2)7^(fc3)) = (2^)' Bnih, k2, ks) S(ki + fc2 + fca), (1.2) 

where (27r)^ is a normalization factor.^ The Dirac delta function imposes a triangle con- 
dition onto the wavevectors fci, ^2, and fcs, ensuring momentum conservation as required 
by translational invariance (homogeneity and isotropy of the universe). As a consequence, 
the bispectrum depends only on the magnitude of the wavevectors, the wavenumbers fci, A;2, 
and ks- To simplify the analysis, we introduce a normalized, dimensionless shape function 
S{ki, k2, ^3), where an overall scaling with k~^ is removed, 

Sih, k2, h) := ^ ikik2hf Bniki, k2, ks). (1.3) 

^Note that Bn is sometimes called F, for example in [49]. 
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(a) tetrapyd domain Vt (b) NG shapes 



Figure 1. (a) The tetrapyd domain Vj- as given by eq. (1.4). (b) Shape templates, taken from [50], 
figure 30. 

Here, N is chosen such that S{k,k,k) = 1; this can be achieved for the most common 
scale-invariant templates with N = const. 

The momentum space domain of S is tetrahedral, but given the existence of a minimally 
observable length scale, the domain is capped by some fcmax resulting in a "tetrapyd" , 

Vr : ki < k2 + ks for ki > k2, ^3, 
A;2 < fci + ks for k2 > fci, ^3, 

(1.4) 

ks < ki + k2 for ks > fci, k2, 
^1,^2, ks < k 

max 5 

see figure 1(a). 

1.2 Shape classes 

Do certain classes of shapes exist? Since we will encounter many shapes in this paper, it 
is prudent to give a brief overview of known shapes, following [30]. In figure 2 we plot 
three illustrative shape functions 5(1,^2/^1,^3/^1). Because of the symmetries of S on 
the tetrapyd domain, it is sufficient to plot momenta satisfying ks/ki < k2/ki < 1, see 
eq. (1.4); the triangle inequality requires ^2/^1 + ^3/^1 > 1. All models discussed here are 
scale-invariant. 

First, there are center- weighted models that peak on equilateral triangles. The most 
prominent shape function in this class is the equilateral template S^equi, see figure 2(a), which is 
an approximation to the dominant bispectrum contribution in DBI inflation [15, 51] (reviewed 
in [52]). We give an explicit formula for the shape in eq. (2.21). The orthogonal shape 
proposed by Smith et al. [49] is roughly given by the equilateral shape with a constant 
removed, rendering it similar to the so called enfolded shape proposed in [53] (see also [54, 55] 
for closely related shapes.). 

A second class encompasses corner-weighted models, which peak on squeezed triangles 
(ki <^ k2 — ks and permutations). They feature a divergence in the ki ^ limit, as seen in 
figure 2(b). This local shape is featured in many multi-field models [56] but also more exotic 
proposals such as the ekpyrotic scenario [57] among others; it can be recovered easily via the 
non-linear 5A^-formalism [13, 14]. 

We complete our brief overview with edge-weighted models that peak on flattened tri- 
angles. There are only few examples known to have a dominant contribution of this type. 
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(c) *S'single 



Figure 2. Plots of some illustrative shapes as introduced in [50]; see figure 1(b) for interpretation, 
(a) The equilateral shape, (b) The local shape that peaks on squeezed triangles, (c) The so-called 
"single" shape S'singie = ^111/^3 that is non-zero for enfolded triangles. 

for example models with non-Bunch Davies initial conditions [58]. They also constitute one 
of the non-suppressed but unobservable small contributions in standard slow-roll single-field 
models. Fergusson and Shellard call this contribution the S^singie template, plotted in fig- 
ure 2(c). We address this shape in section 4.2.1. 

Various other shapes exist in the literature, some of which include explicit scale de- 
pendence or localized features, highlighting the need for a unified approach to accommodate 
general shapes. 

2 Separable, polynomial expansion of the bispectrum 

To deal with general primordial shapes,^ we briefly review the separable mode expansion of 
the bispectrum using tetrahedral polynomials, as introduced in [27, 29]. Our goal is to analyze 
a general shape function S{ki, k2, ks) by a decomposition into orthonormal or separable mode 

^We do not perform a decomposition in spherical harmonic space, which requires an evolution of the shape 
to later times; the latter is needed to make direct contact with CMB observations [29]. We give an outlook 
to the late-time CMB bispectrum in section D. 
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functions. The definition of either basis is explained in section 2.1. The conversion matrix 
between respective expansion coefficients is provided in eq. (2.20). 

Since we need to compute the correlation between shapes, we start be defining an 
integration over the tetrapyd domain Vj- 



T[f] := / /(fci, A:2, h) w{ki, A:2, ^3) dVr, (2.1) 
JVr 



employing a weight function w{ki^ A:2, ^3) that should be chosen to reflect the scaling of the 
CMB bispectrum estimator. To remain compatible with work done by Fergusson, Shellard 
and Meerburg, we employ w = 1 for the analysis on the primordial wavenumber domain.^ 
An inner product is imposed by letting 

{f,9):^r[f9% (2.2) 

where a star denotes complex conjugation (not needed for polynomials, but required for the 
Fourier modes of [37], see section 5.1). The correlation between two shapes and 5^2 can 
then be measured by the normalized inner product on the tetrapyd 

:= , ^^'' ^ '^ (2.3) 

In [59] and related works, this correlation coefficient is called the "cosine" between shapes. 

2.1 Construction of an orthonormal basis 

Let us rescale wavenumbers by^ 

and generate orthonormal polynomials in a one-dimensional subspace as a starting point. To 
this end, we define the reduced weight function 



w{x) w{x,X2,X3)dVr,x2,x3i (2-5) 

where we integrated out X2 and X3, as well as a reduced, one-dimensional integral 

%[f] t f{x)w{x)dx. (2.6) 
Jo 

Using w{x^ z) — 1, we get w{x) — ^x(4 — 3x). As the basis for the polynomials we choose 
monomials, 

Cr := r G Nq. (2.7) 



^For the CMB bispectrum, the weight in [27], eq. 48, should be employed. Note that there is a mistake in 
one of the equations in [27] that we correct in appendix D.l. In order to prepare the analysis of primoridial 
bispectra for the multipole expansion, one should use the weight w = l/(A:i + /c2 + fe) as mentioned by 
Meerburg in footnote 4 of [37] . 

^Unless scale-dependent bispectra are discussed, the maximal wavenumber is irrelevant. 
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Figure 3. Plot of the first ten one-dimensional basis functions, qr{x) (compare to [27], figure 5). 



Defining 



Wr 



r + 6 



2(r + 3)(r + 2)' 

we arrive at orthonormal polynomials (w.r.t. the reduced weight function) via 



(2.8) 



qr{x) := 



Wo 


Wl . 


Wr 


Wl 


W2 . 




Wr-1 


Wr . 






ei . 





(2.9) 



We choose J\f such that 7^[g^] = 1 for all r, so that the polynomials qr{x) fulfill 

{qr,qs)x TxiQrqt] = ^rs' 



(2.10) 



Analytic expressions for the first five functions can be found in [27], eq. 54. See figure 3 for 
a plot of the first ten orthonormal polynomials. 

Due to the permutation symmetry on the tetrypyd domain, we define a three-dimen- 
sional complete set of functions as a symmetrized superposition 

Qn{x, y, z) := q{pqrqs} = ^ [qp{x)qr{y)qs{z) + qr{x)qs{y)qp{z) + qs{x)qp{y)qr{z) 

+ qp{x)qs{y)qr{z) + qs{x)qr{y)qp{z) + qr{x)qp{y)qs{z)] . 

The definition of the super-index, the mode number n, is arbitrary;^ we map n to the triplet 
{p, r, s} using the "slicing" ordering, where the {p, r, s} are sorted in ascending order of the 
quantity + + r^. To clarify the ordering, we print the n ^ {p, r, s) mapping up to n = 53 
in appendix A. We choose Af such that T[Qr] = 1 for all n. 

The number djy of independent symmetrized products Qn{x,y,z) of polynomial order 

N is 



{(In} = {1, 1, 2, 3, 4, 5, . . .}, c^at = 1 + dN-2 + (Ins - (Ins- 



(2.12) 



Thus, we have 53 mode functions if we include all polynomials up to order 9. These polyno- 
mials are separable, but not orthonormal; nevertheless, their inner products 



In 



(2.13) 



^One may want to investigate a different ordering to improve rapid convergence for common shapes. The 
chosen order has proven sufficient, but is not the result of any optimization. 
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Figure 4. (a) Plot of the inner product matrix F defined in eq. (2.13) of the first 31 modes (compare 
[27], figure 7). (b) Plot of the conversion matrix A' defined in eq. (2.15) for the first 31 modes 
(compare [27], figure 7). 



show a predominantly diagonal structure, inherited from the orthonormal one-dimensional 
functions. In figure 4(a) we plot the color-coded matrix components for the first 31 mode 
functions (up to polynomial order 7). The "checker-board" pattern stems from sets of poly- 
nomials of the same order with lengths in eq. (2.12). It is worth noting that the Qn become 
more correlated the higher the mode number n is. 

In a third step we obtain an orthonormal basis {Rn}^ 

{Rn,Rm) ^nm, (2-14) 

using Gram-Schmidt orthonormalization. Let 

-{Qn.R'm)^ n>m 



nm 



1, n = m , (2.15) 

0, otherwise 



where 

n 

< = E ^nmQm. (2.16) 
m=l 

We obtain the final orthonormal functions by letting 

Rn := , (2.17) 

The color-coded matrix entries of A' for the first 31 mode functions can be found in figure 4(b). 
2.2 Mode expansion of bispectra 

Equipped with the symmetrized, orthonormal mode functions on the tetrapyd, we can express 
the bispectrum Bn of any theoretical model by a series of mode coefficients an of the shape 
function 5(^1,^2,^3) in eq. (1.3). The full shape function is recovered as the formal sum 

00 

S{ki,k2, ks) = ^ Rn{xi,X2, X3), (2.18) 

n=l 
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Figure 5. Plot of the expansion coefficients a^^"^""'^ and a^^^^^ equilateral shape (blue), as 

defined in eq. (2.21), and the DBI model (green), as defined in eq. (2.22); compare [27], figure 10. 



where the momenta are rescaled as in eq. (2.4) and the coeflficients are given by 

a^ = {Rn,S). (2.19) 

The matrix 

A=(A.^), Xr,m= /7^^ ^ (2-20) 

can be used to perform the conversion between the expansion coefficients and the coeffi- 
cients of the expansion in the separable Qn basis, denoted a^. Since the conversion matrix 
needs to be computed only once, we retain the full advantage of separability in computa- 
tions, given that the expansion series is sufficiently convergent (see section 2.2.1); in this 
manner, shape functions can be expressed by a few numbers. In the following, we obtain 
mode coefficients numerically up to rimax = 53 for several shapes.^ 

To illustrate the modal decomposition, let us investigate the shape originating in DBI 
inflation [15, 51]. This stringy inflationary model incorporates a speed limit in fields space 
via a modified kinetic term in the Lagrangian and is a special case of the more general 
Lagrangians that we discuss in section 4. The resulting bispectrum led to the definition of 
the equilateral template. Thus, the expansion coefficients of these two shapes, namely^ 

5equi(fei, k2, h) := ^ (i^i2 - - 2Kni) , and (2.21) 
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5„e,fe.fe,«.= ^(-f + §| + f), (2^22) 

should have a strong correlation. The shorthand notation A^,,, indicates simple polynomials 
in {ki} and is explained in appendix B. We find the mode coefficients shown in figure 5, re- 
covering the analysis by Shellard et al. How good, quantitatively speaking, is the convergence 
and correlation between these shapes? 



^We only use the coefficients in the text and figures. Hence, we drop the superscript 7^ on a in favor 
of a superscript to indicate the shape function used. 
^Equivalent expressions used in the literature are 

n (K — 2k ) 1 

>Sequi(A:i, A:2, ks) := and 5dbi(A:i, /cs, fe) = ^ (i^s + 2Ki4 - 3K23 + i^iis - ^Ki22) • 

Km AiiiA^ 
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Figure 6. Plot of the sum of squared coefficients in eq. (2.25) for the equilateral shape and the DBI 
model (see [27], figure 9). The inset shows that O (15) modes are sufficient to achieve better than 
99.5% convergence. 



2.2.1 Convergence properties 

We would like to use the partial sum 



'Sn.ax E (2.23) 



n=l 



to approximate a given bispectrum. In order to quantify the level of convergence of that 
sum, we use Parseval's theorem, 

oo 

= (2.24) 

n=l 

and investigate how fast the limit 

anC^n — ^ 1 (2.25) 



is approached. 

We consider a series to be sufficiently convergent if the series truncated at order rimax 
reaches 1 — ^, where 5 is a small parameter. Our results for the mode expansions of the 
equilateral shape and the DBI model are shown in figure 6. Using rimax = 53 modes, we get a 
convergence as good as £ ^ 0.0002. Fifteen modes are sufficient to achieve e ^ 0.0031. In the 
literature, a minimal convergence of 95% to 98% is required [29], but a higher convergence 
may be needed for comparison with PLANCK data. 

2.2.2 Correlation between shapes 

The correlator between two shapes is defined in eq. (2.3). For the equilateral shape and the 
DBI model we get Cequi,DBi ^ 0.991391, an expectedly high correlation. In figure 7, we show 
the partial correlation sum 

n 

^^(e<iui)^(DBl)*^ n<n_, (2.26) 

i=l 

which converges to the final result Cequi,DBi rather quickly, because the sum over the squared 
mode coefficients in eq. (2.25) converges satisfactorily after only 15 modes for both shapes. 
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Figure 7. Plot of partial correlation sum for the equilateral shape and the DBI model. 

We observe that certain modes improve correlation more than others, which is expected 
since no physical input went into the ordering of the basis functions. If we had a particular 
shape to match, we could reorder modes to achieve faster convergence. Such a procedure 
should be performed if observational data (CMB or LSS) is analyzed, in contrast to the 
approach in [29] which used the same polynomials as above. 

2.2.3 Classes of shapes by correlation 

Given the polynomial decomposition, which other shapes have good convergence properties? 
Fergusson and Shellard [28] give a selective overview of common bispectrum classes. We 
find good convergence for the center-weighted models, such as the DBI shape, that peak on 
equilateral triangles; this is not surprising, since all models in this class are highly correlated 
with the equilateral shape Sgqui- 

The expansion series of corner- weighted models that peak on squeezed triangles do not 
converge due to the divergence of these models in the ki ^ limit. Fergusson and Shellard 
[28] suggest a clipping of such shapes based on a minimal physically relevant wavenumber 
kmin ^ (2//max)fcmax- Commonly discussed shapes in this group are correlated among each 
other by varying degree. It should be noted that the local shape is also correlated with the 
equilateral shape and not orthogonal to it. Thus members of this class are often distinguish- 
able, but not completely independent. 

The correlation of shapes among members of these classes is summarized in table 1, 
whose entries stem from table I &; II in [30] and were reproduced using our framework. 



Model 


equilateral type 


local type 


feature type 


equilateral type 
local type 


at least 86% 


up to 46% 
up to 62% 


up to 36% 
up to 44% 



Table 1. Shape correlation according to eq. (2.3) among members of three shape classes: the equi- 
lateral template, the DBI model, and the ghost model make up the class of equilateral-type NG. The 
class of the local-type NG consists of the local template, the warm model, and the flat model. The 
feature model is defined in eq. (2.34). Values are take from [30], tables I & II, which we reproduced. 
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Figure 8. Plot of the warm shape, where the cut-off (top left corner) in eq. (2.31) is marked by a 
red line. 

2.3 Divergent shapes and the clipping mechanism 

The shape derived from the local ansatz,^ 

1 Ks 

S'locai o — ' (2.27) 
6 Km 

does not have a convergent expansion series due to the divergence seen in figure 2(b). To ame- 
liorate similar problems with divergent shapes, one can try to remove any local contribution 
from the discussed shapes or, if a divergence remains, clip the shape. 

Several corner- weighted shapes that peak for squeezed triangles are discussed in [28], 
section III.E. For example, the warm shape [60] produced during warm inflation (inflation in 
the presence of a radiation bath, see [61] for a review), 

S'warm OC — , (2.28) 

^333 

is one of the shapes that require a cut-off in order to be treatable by the mode expansion 
algorithm. Based on the flat sky approximation 

^ , (2.29) 

^0 - m 

where tq is the age of the universe and tr the time from recombination till today, Fergusson 
and Shellard [28] argue that a natural cut-off exist at 

2 

fcmax- (2.30) 



Accordingly, the divergent part can be removed by requiring [28] i.e. 



k2 + k'i 



< 0.015, (2.31) 



where the exact value of the cut-off is chosen by hand. The contour of this cut is marked red 
in figure 8. We provide a discussion of the clipping mechanism's validity in appendix C. 



^The local ansatz is defined in eq. (3.19). 

^°Next to the warm shape, non-Gaussianity of the local/constant type is also produced if damping terms 



are allowed to depend on the temperature [62]. 
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Figure 9. (a) 3D plot of the feature model's shape function for = k^a^/ojf^ where cjj = 9, and 
(5 = 0. (b) Plot of the sum of squared coefficients for the feature model using a polynomial expansion 
(see [37], figure 3); convergence is slow (compare to figure 27, where faster convergence is achieved by 
using Fourier modes). 



However, the warm shape is still hard to treat with the modal expansion. First, the 
effects of the cut-off need to be kept small. This can be achieved, as proposed in [28], by 
smoothing the shape after truncation, i.e. by applying a Gaussian window function. Second, 
the divergence of the warm model in the limit A:3 ^ has a direction-depending sign, contrary 
to the local shape in eq. (2.27), to which it is only 33% correlated. 

In the following, we will encounter shapes that are likewise divergent, but most of them 
are less problematic: subtraction of the local contribution, 

^(conv) _ ^(div) _ sf^l^^Slocah (2.32) 

renders these shapes convergent. Here, we defined 

(2.33) 



/nl 



k2=ki 



which is an analog of the non-linearity parameter fj^i^^^ for shape functions S}^ 
2.4 Scale-dependent features 

Bispectra with scale-dependent features occur occasionally, for example the sharp feature 
models 

S'^^\h, k2, h) = Sin (^^^1^ + 5) , (2.34) 

found by Chen [32, 33], see also [63]. Here, A:* defines the scale of horizon exit when the 
inflaton field encounters a feature in the potential, and 5 is a phase. 

The convergence of such oscillatory bispectra w.r.t. the polynomial expansion is rather 
slow, see figure 9(b). Meerburg [37] proposed an alternative mode expansion using Fourier 
modes for such oscillatory models, to which we come back in section 5. He found faster 
convergence of the expansion series for certain models such as the feature model, see figure 27, 
resonant type non-Gaussianities [32, 64, 65] and models with non-standard vacua [53, 66]. 



^^Note that /nl^ and /n£^ differ by a shape-dependent multiplicative factor due to the normalization 
S{k,k,k) = 1. 
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These shapes highhght a common problem of truncated modal expansions: the choice 
of basis functions renders a truncated expansion more sensitive to certain shapes, while 
becoming blind to others. We investigate an example of a problematic shape in section 5.2 
that cannot be recovered rapidly with a polynomial expansion, while Fourier modes work 
well. 

3 Review of general single-field infiationary models with non-standard ki- 
netic terms 

In this section we review the analysis of general single-field models based on the most general 
Lagrangian that is Lorentz invariant and whose kinetic term contains first derivatives in the 
field only. The reader might want to skip sections 3.1 and 3.2 if he/she is familiar with 
general slow-roll single-field models. In addition, readers familiar with the slow-roll violating 
models by Noller and Magueijo [23] might want to skip section 3.3. 

The single-field inflation model currently favoured by data, i.e. the simplest model 
that is not excluded by observations, has a standard kinetic term, is minimally coupled to 
gravity and follows a slow-roll trajectory. Such a model only produces a negligible amount 
of non-Gaussianity [18]. Since future observations [1] may detect a NG signal (note the 
nearly 2a signal in eq. (1.1) for a plethora of more complicated inflationary models 

allowing for larger non-Gaussianities have gained popularity. For example, the DBI model, 
as introduced above, invokes a non-standard kinetic term and naturally leads to observable 
non-Gaussianities. A general Lagrangian encompassing DBI inflation among other models 
was proposed by Garriaga and Mukhanov in the context of k-Inflation [67, 68]. Here, a 
speed of sound different from unity is possible so that slow-roll conditions of the field can be 
violated without spoiling the flatness of the power spectrum. These are the models we wish 
to discuss in section 3.3. 

Even more general models, dropping the requirement of first order field derivatives in 
the Lagrangian in favour of at most second order equations of motion [69], go under the 
name of Horndeski models. Of particular interest are Galileon models [24, 70] where the 
self-interactions are protected by a covariant generalization of the Galilean shift symmetry 
(corresponding to a specific choice of coefficients in the Horndeski action); some of these 
were investigated by Seery and Ribeiro in [24], who found a "new" Galileon shape which 
is however strongly correlated with a known shape appearing in radiatively stable Galileon 
models with higher order derivative operators (discussed in [54]), and to a lesser degree to 
the one found in [55] brought forth by second order slow-roll corrections. 

We focus on the narrower class explained above, with the aim of providing a complete 
anatomy of the resulting bispectrum. To this end, we slowly increase the degree of complexity, 
introduce key concepts and set our notation in three steps: first, we summarize the framework 
of Seery and Lidsey [21] in section 3.1, who calculated the bispectrum for a general single- 
field scenario at first order in slow roll under the assumption that the speed of sound is close 
to one. They found that the bispectrum can be expressed in terms of the standard slow-roll 
parameters and three additional parameters that measure the non-trivial kinetic structure, 
the speed of sound and the rate of change of the speed of sound. 

We follow with a summary of Chen et al. [16, 22] in section 3.2, who treat a more 
general case, sometimes referred to as slow variation in contrast to the conventional slow-roll 
case. They calculate the bispectrum to leading order in slow roll, but allow for Cg <C 1, 
and find a bispectrum depending on five free parameters. The interesting case C5 <C 1 can 
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generate non-negligible amounts of non-Gaussianities, leading to a bispectrum dominated by 
two possible qualitative shapes; the corresponding amplitudes depend on two parameters. 

Lastly, we focus on the work by Noller and Magueijo [23] in section 3.3, which is moti- 
vated by the insight that deviations from the usual slow-roll conditions do not always break 
scale-invariance of the power spectrum [71, 72]. Their maximally slow-roll violating models 
allow up to e ^ 0.3 consistent with observations. 

The violation of slow roll has also been investigated in the case of the Galileon models 
{G-inflation, [73]). 

3.1 Non-standard kinetic terms in single-field infiation assuming slow roll 

Following [21] we treat a single- field inflationary model that allows for non-standard kinetic 
terms by letting the Lagrangian £ be a general function of the form 

C^P{X,<l)), with X = -/^V„</.Vfc</>, (3.1) 

SO that it contains at most first derivatives in the scalar field and is Lorentz invariant. The 
field (f) is minimally coupled to gravity leading to the action 



The scalar field has a speed of sound 



2+^ 



(3.2) 



which can differ from one. Here, (. . . denotes a partial derivative w.r.t. X. 

3.1.1 The slow-roll approximation 

Slow roll can be imposed by requiring^^ 

These flow 'parameters describe the evolution in parameter space of inflationary models. We 
decompose the first flow parameter as 

4> OH X dH 

'^'iP^'lPdx^'-'^^'''- ^^-^^ 

Smallness of e^^ induces 0^ <C . With these definitions, the field equation leads to X = 
— 6i7c^X(l — e^/e), enabling the replacement of derivatives w.r.t. t with derivatives w.r.t. 
the field's kinetic term X. 

The speed of sound is best described in terms of an additional parameter 

u:^l-c-^^-2X^, ^^-4-- (3-6) 



^^Often an alternative definition r/H — 2H" / H is used, which corresponds to riv — V" /V for single-field 
models with a canonical kinetic term, leading to 77 = —2riv + 4ey at first order in slow-roll. 
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Here, we introduce a new "slow-roll" parameter s to express the rate of change of C5, so that 
u = 2Hs{l — u). We require the speed of sound to be slowly changing in order to retain an 
almost scale-invariant spectrum. In [21] e and 77 are assumed to be of the same order (e ^ 77) 
and constant, consistent with e,r] ^ O (e^). Demanding u = O (e) implies a small deviation 
from = 1 leading to 5 = O (e^); it follows that ex = O (e^) [21] so that ex is subdominant 
to e(f). To simplify the description of the bispectra in this framework, we define 

E:^XPx + 2X^Pxx, (3.7) 
A := X^P^xx + Ix^'Pxxx , (3.8) 



which can be rewritten as 

" H\{l-u) • ^ 



- — {l-u)s-u 
6 ex 



c1 " " "6 

These identities hold in general and are not contingent on the slow-roll approximation. 



(3.9) 



3.1.2 The power spectrum 

Seery and Lidsey [21] analyze the action S[7l] for the comoving curvature perturbation TZ to 
leading order in slow roll using the ADM formalism. Keeping only terms up to quadratic 
order in TZ in the action, the Fourier transform of the two-point correlation function, the 
power spectrum, reads 

(7^(fcl)7^(fc2)) = {27Tf6{ki + k2)P{ki), (3.10) 

m = fp. (3.11) 

The primordial power spectrum becomes 

V^(k) := ^,P(k) = (3.12) 

and the scalar spectral index can be computed to 

dlog'P7^ ^ o o ^Q1Q^ 

ris-l —- — 2^ -2e - 7/ = -2e^ - 2ex - V- (3.13) 

d log k 

3.1.3 The bispectrum 

In [21] the bispectrum is calculated by expanding the action up to third order in 7^, restricted 
to leading order in slow roll. Seery and Lidsey apply a path integral formalism for the three- 
point function, 

/[d^] exp [iJcL) 

and analyze the tree level contributions from all possible interaction vertices (7^^, TZd~^TZ, 
lZd^TZd~^lZ) . Combining these terms, the Fourier transform of the three-point correlation 
function is (assuming ki ^ k2 ^ ks) 

{n{k,)n{k^)n{k,)) = (27r)^ 5^^) (E, fc,) rl^^,, (3.15) 
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where ([21], eq. 82f) 

A = —{u + e)K22 - ;p ( ^ + ^ 3 ) ^222 - ^^23 + "(t/ - 1^ - 6)^3 + -^^12- (3.16) 

We employ the notation introduced in appendix B, which simphfies the identification of 
shape templates in the bispectrum amplitude. We can construct a shape function as defined 
in section 2 by letting S — A/ {Km N)^ where N should be set such that S{k^ k, k) = 1. 

As an aside, we would like to comment on the validity of the consistency relation found 
by Creminelli and Zaldarriaga [74]. In the squeezed limit ks 0, which corresponds to a 
renormalization of the background, 

(7^(fcl)7^(fc2)7^(fc3)) = i^Trf 5^^^ ki) (1 - ns)P{ki)P{h) (3.17) 

is valid for general Lagrangians [75], and initial conditions [76]. This consistency relation 
shows the existence of local-type NG for any inflationary model with 7^ 1. The main 
reason for the existence of this relation is the presence of a single dynamical degree of freedom 
that acts as a physical clock for exiting the quasi-de Sitter phase. 
The amplitude A in the e = 77 = case, namely eq. 88 in [21], 

^U.. = 4„§J-4„§|2-2„:^-|a'3, (3.18) 

highlights the NG contributions that arise purely from the inflaton field, i.e. without the 
coupling to gravity. In particular, A\^^^^q ^-^^ follows in line with the consistency relation 
(independent of the speed of sound). 

3.1.4 Bispectrum amplitudes and /nl 

To search for the presence of non-Gaussianity, it is useful to define a single quantity which 
provides a measure of the overall amplitude. Naturally, any shape and k-dependence is 
lost by focusing on just one number (see [77] for a proposal to go beyond the amplitude by 
considering the running of the bispectrum). The most common approach consists of picking a 
separable shape template, such as the local or equilateral one, whose magnitude can be easily 
extracted from the data. Theoretical bispectra are then projected onto the template to enable 
a comparison with the measured value, which can be problematic if the correlation is weak. 
A better definition of an amplitude that does not rely on templates was recently employed in 
[29] and uses the modal decomposition in spherical harmonic space to compare theory with 
observations. However, since the latter approach is still being optimized (choice/ordering 
of basis functions) and the forward evolution of shapes by means of the radiation transfer 
functions is non-trivial, we shall use the standard definitions of the non-linearity parameters, 
as defined in eqs. (3.21) and (3.22), at this point. 

Consider first the so-called local non-Gaussianity, corresponding to fci ^ A:2 ^ fcs, 
whose amplitude is expressed by the non-linearity parameter f^if^. It was introduced in a 
(spatially) local ansatz for non-Gaussianities [10-12] caused by the square of the Gaussian 
scalar potential^^ 

n^nG + p^L{nl-{nl)). (3.19) 



^^Note that the sign convention of Seery and Lidsey [21] is opposite to the more common one by Komatsu 
et al. [19]. We employ the latter in order to keep comparability with the WMAP results. 
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The corresponding three-point correlation function is 



(7^(kl)7^(k2)7^(k3)) = (2^)^ 5^ j^/nl ^^(fc)^- (3-20) 

Given the general single- field bispectrum to leading order in slow roll in eq. (3.16), the 
prediction 



rlocal ._ 5 A 



12 



-^(ris-1) (3.21) 



is a direct consequence of eq. (3.17). 

The other commonly used limit is the equilateral one fci ^ /c2 ^ fca, leading to the 
prediction [21] 



requi 



A 



6 E.fcf 



275 10 e 55 5 

u s + + ^77. (3.22) 

k.=k,=ks 972 729 ex 36 12 ^ ^ ^ 



This reduces to Maldacenas [18] result f^' = -^(n^ - 1) - ||(n^ - 1) for = 1. 

Thus, the overall amplitude of non-Gaussianities for these templates is of the order of 
the slow-roll parameters and not observable, because foreground effects are expected to yield 
an |/nl| of C^(l) to (9(10) (see [78, 79]). In reaching this conclusion, we assumed that at 
least one (or both) of the templates have a good correlation with the primordial shape, which 
is indeed the case as we shall see later on. 

3.2 Towards a violation of slow roll: a small speed of sound 

The discussion in the previous section on general single-field, slow-roll inflationary models 
with cl ^ 1 showed that non-Gaussian signals are small, because f^^ is slow-roll suppressed. 
But we have also seen in eq. (3.18) that a pure field contribution arises for < 1. Further, 
the DBI model in section 2.2 gave a first indication that a low speed of sound can indeed 
lead to observably large non-Gaussianities. Thus, let us turn our attention to models with a 
small speed of sound, following [16, 22]. 

3.2.1 The power spectrum 

In the more general class of models in [16], the power spectrum has an additional scaling 
oc in comparison to the result of Lidsey and Seery in eq. (3.11),^^ 

Here, the (sound) horizon exit takes place at Cgk = aH, and the spectral index 

ris-l^ -2e-r]-s (3.24) 

has an additional dependence on the variation of Cg via s. The tensor-to-scalar ratio is also 
altered to r = Vt/Vc^ — —ScsTIt = 160^6, where ut is the tensor spectral index. Evidently, 
a deviation of the speed of sound from unity enhances the amplitude of the scalar curvature 
fluctuations relative to the amplitude of tensor fluctuations. 



^^In order to keep our equations comparable to the original work, we express the spectra in terms of 
the curvature perturbation on uniform density hypersurfaces ( in this section; the latter is identical to the 
comoving curvature perturbation 71 on super horizon scales. 
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3.2.2 The bispectrum 

In order to keep the power spectrum almost scale- invariant, Chen et al. assume that all 
of the slow-variation parameters 6,77,5 and / in eq. (3.32) are of the same order O {e) and 
small e <C 1. Using Maldacena's approach as in [21], they derive the three-point correlation 
function from the effective action of the curvature perturbation ( expanded up to O (C^)- 
After a field redefinition to absorb a term proportional to the variation of the quadratic 
Lagrangian, they find that for an arbitrary speed of sound some leading order terms are not 
slow-roll suppressed and, thus, potentially observable. Since these terms vanish in the limit 
C5 ^ 1, subleading corrections need to be kept. By these means, they obtain the complete 
bispectrum up to O (e) 

(C(ki)C(k2)C(k3)) = {2nyS^{ki + k2 + ks)Pl{K)^ 

X {Ax + Ac + Ao + Ae + Arj + As) , (3.25) 
where the bispectrum amphtude has been decomposed into 

A=(i-l-|[2-(3-2c,)q)^^. (3^26) 

(3.28) 



A = .f-^ + ^ + ^y (3.29) 



K , 

A,^V[^], (3.30) 



As = sFs. (3.31) 

Here, Ci is the Euler constant. The shape functions F.. in Ao and As are given in appendix 
B.l of [16] as unevaluated integral expression; exphcit formulae can be found in appendix D 
of [70] . The definitions of the four slow- variation parameters are 

'-■^-IP^ ^^^A^- ^^-^'^ 

The additional parameters A and S are defined in eqs. (3.7) and (3.8). It is worth noting that 
the amplitudes of the bispectrum shapes are determined by 5 parameters only: e, 77, 5, Cg and 
X/Tj. The parameter / is only contained in subleading contributions (see [16], section 4.4) 
leading to O (e^) effects. 

Pc,{K) is the power spectrum in eq. (3.23). Note that the definition of the bispectrum 
amplitude A by Chen et al. differs by a factor of 4 from the one by Seery and Lidsey in [21]. 
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3.2.3 Bispectrum amplitudes and /nl 

Chen et al. [16] use the same definition for non-linearity parameters as in eq. (3.20); focusing 
on equilateral triangles (ki ^ k2 ^ ks) gives-*^^ 

^ + (3-2c0-^ 

(3.33) 



/nl 


81 U 


- 1 - 


/nl 


35 / 
~ ~108 \ 




/nl 






unl I 







We can learn from these amplitudes that observable non-Gaussianities can arise in these infla- 
tionary models if Cg is much smaller than one, and/or A/S gives a non-negligible contribution 
to /nl- 

Can we already constrain Cg based on current limits on non-Gaussianities? Senatore 
et al. [49], who employ a slightly different formalism for their general single-field inflation 
model, use the equilateral and an orthogonal shape orthogonal to the equilateral shape w.r.t. 
to the scalar product in eq. (2.3)). Using these separable shapes, they were able to provide 
(weak) constraints on the speed of sound: > 0.011 at 95% CL or is so small that the 
higher-derivative kinetic term dominates at horizon crossing. 

3.2.4 Recovering slow-roll results 

In the case of conventional slow-roll inflation, i.e. for <C 1, Chen et al. recover the result 
-4sl/4 (the factor of 4 is due to different conventions) from eq. (3.16). To this end, they 
neglect As and Ao^ which are second order in the conventional slow-roll parameters. 

+ e - — + -r— + —r;- ] + V 



Furthermore, for u = and 5 = (i.e. = 1 = const.), they find agreement with Maldacena's 
result up to another prefactor of 8:^^ 



A\ A 3nr-2n, -1 1 - ny 



(3.35) 



16 16 
3.3 Allowing as much slow roll violation as possible 

Noller and Magueijo took inspiration from work by Khoury and Piazza [72], who derived the 
bispectra for exactly scale-invariant spectra without restricting the parameter space by slow 
roll. Dropping the requirement of exact scale invariance, Noller and Magueijo [23] compute 
the general single-field bispetrum using the variables 

-,26 + 65 26 - es . . 

ai:= 71^-1 = -, a2 := ^ r, (3.36) 

6. + 6-1 6. + 6-1 



^^Note that there is a mistaken parenthesis in [16], eq. 5.3. Note that the sign must be inverted for 
comparabiUty with the original work. 

^^Due to yet another convention, there is an additional factor of 2, so that in total, Maldacena's convention 
and the one used here differ by a factor 8 
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where = 5 = Cs/{csH) and instead of ex they employ 



fx:-- 



3ex' 



(3.37) 



Note that the parameter A can be written as 



(3.38) 



according to eq. (3.9). They assumed that e and dominate over 77 and 77^ = es/(esH) (as 
weh as higher order parameters), which are set to zero in the fohowing. 

3.3.1 The bispectrum 

Fohowing the same hne of thought as outhned above, an effective action of cubic order in 
can be obtained, but in contrast to Chen et al., no slow- variation conditions are assumed. 
Computing the full bispectrum in this way for each term in the action, they obtain 



A, 



1 



4c2 V 2K^ J 



Us — I 



Us — I 



(€ + €,- - l)X.3(a2) + c2(e + e, - l)X.3(ai) 



(e - 3)X^^2(q:2) + Sc^X^^zlai 



ria — l 



A 



1 



Us — I 



Us — I 



A,2 = 



16c2 



(3.39) 



where 



%(o;) = cos ^r(3 + 



2^^2(0:) = cos — r(l + a) 



(2 + a 



(1 + a 



2^C(9C)2 (") = - cos ^r(l + (x)K-2 



K 

K Kii 



K23 



+ K: 



111- 



1 + a 



COS ^r(l + a) 



^CeC^x^"^ ^ — r(l + a 



a - 1 K 

1 — a A A ^ 1 — a 

+ ^—^12 - 2-^i^23 - 2aKui 



K12 - cxKin 



Z,2(a) = cos — r(l + a)(2 + a/2) [K3 - i^i2 + 2Kiii] . 
These amplitudes reduce to the slow-roll limit discussed in the previous sections. 



(3.40) 
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3.3.2 Bispectrum amplitudes and /nl 

The non-linearity parameter for the equilateral shape is given by 

A 



/nl = 30 ^3 



(3.41) 

ki=k2=k2, 



equivalent to the definition in eq. (3.22) (up to a factor of 4, because Noller and Mageuijo's 
amplitude A corresponds to Chen et al.'s convention). The full result can be expressed as 
[23] 

/nl = Ci(n„ e) (1 - C2(n„ e, fx)c-^) , (3.42) 

where the Ci are functions of the slow-roll parameters and fx (given in the appendix of [23]). 
In [23], the dependence of /nl on the parameters is discussed for some examplary models. 
For instance, the case /nl ^ 1 leads to 

/nl~0(c;2)+o(^). (3.43) 

The running of the non-linearity parameter 

_ dlog|/NL| , . 

dlogfc ' ^^-^^^ 
derived from the parametrization in eq. (3.42), is investigated in detail in [23]. 

4 Correlation properties and mode decomposition of bispectrum shapes 
in general single-field models 

This section contains our first main results. We analyze the basic terms that arise in the 
bispectrum amplitudes of single-field models, followed by a discussion of some concrete fast- 
roll models. 

4.1 Parameter extraction and the CMB bispectrum 

The motivation behind a series expansion of the bispectrum is the easier comparison to the 
data for general not necessarily separable shapes. To this end, the predicted CMB bispectrum 
Bi^l^l^ needs to be expanded in terms of the multipole moments (not wavenumbers) for each 
theoretical model (after evolving it forward in time using radiation transfer functions) and 
compare the resulting coefficients to the observed ones. Based on WMAP5 data and the 
polynomial modal decomposition, Fergusson et al. [29] were able to constrain the amplitudes 
for certain commonly discussed models and plan to incorporate the WMAP7 data release 
(and subsequently PLANCK data). It is desirable to optimize the choice of mode functions 
as well as their ordering, i.e. via a principal component analysis, which has not been done yet. 
Therefore, we work with the primordial shapes 5(^1,^2,^3) and discuss if the contribution 
of different shapes could in principle be disentangled and whether a parameter extraction 
appears feasible. 

As a first case study, we compare the general slow-variation single-field shape con- 
stituents in eqs. (3.26) to (3.31). If the series coefficients of an observed primordial bispec- 
trum were available with a high signal-to-noise ratio and the expansion coefficients of the 
relevant shape functions were at hand, the model parameters e, 77, c^, 5 and A/H can be 
extracted by a non-linear regression or using MCMC analysis. 
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4.2 Correlation of reduced, non-divergent shape contributions 

The simplest constituents of the shape functions A that we encountered so far are of the 
form 

wherep > l,g > max{2,2?}; and 1 where t>s>r>l. (4.1) 

Ah of these contributions to the Seery and Lidsey result in eq. (3.16) are divergent on their 
own, except for .Aconst Km and v4singie ^ K222/K^- At this point, bear in mind that the 
conversion rule to the dimensionless shape function gives ^single ^ -4singie/^iii — Km/K^. 
In the following, we switch between A and S depending on which one is more convenient to 
use. 

Before analyzing these basic constituents, we would like to remind the reader that the 
general single-field result does not contain the separable equilateral shape ^equi — ^12 — 
Ks — 2^111, see figure 2(a). It closely approximates, however, the divergent-free part of the 
Ae shape, namely 

Ae - 2K3 oc K12 - 3^3 + 8^. (4.2) 
The previously mentioned DBI model is 

^DBi (X A (X K3 - 8^ + 4^, (4.3) 

which is closely approximated by the equilateral shape (see figure 7). These three shapes are 
examples of linear combinations of the K terms where the divergent contributions cancel. 

4.2.1 Correlation matrix of basic contributions to the NoUer and Magueijo 
model 

The K... terms in eq. (4.1) are not independent of each other, see appendix B.l, but we 
can identify five independent, basic, convergent constituents of the general (scale-invariant) 
bispectrum shapes by Noller and Magueijo in eq. (3.40) in addition to the local (oc K^) and 
constant (oc Km) ones. We choose 

K22 

Al K12, A2 ^const ^111, A3 



^ (4.4) 

— ^2 ' — ^3 , — wHsmgle — ^3 , 

and the local shape as a complete set of such basic terms (up to polynomial order 6), i.e. 
we can express all other terms with the right scaling and symmetry properties as linear 
superpositions of the terms in eq. (4.4). 

Let N [•] he the operator that removes the divergent (local) part proportional to K3, 
subtracts the constant mode from a certain shape and normalizes the result, that is 

f-l-'{N[Ai]) = 0, {A2,N[Ai]) = 0, {N[Ai],N[Ai]) = l, (4.5) 

where the non-linearity parameter of local-type non-Gaussianity is given by 

/nTV) = (-5/6) lim [A{k, k, fc3)/(2fc3 + kl)] . (4.6) 
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(a) N[Ki2\ (b) iV[i^iii] (c) N[K22/K] 




(d) N[K2:,/K^] (e) N[Kq/K^] (f) N[K222/K^] 



Figure 10. Plots of the basic normalized and divergent-free contributions N [Ai]^ where i = 1, . . . , 6, 
which are the result of the N [•] operator, cf. eq. (4.5), on the terms in eq. (4.4). 

Since most of the correlation between v4equi and w4dbi is due to the constant mode, we 
achieve better shape discrimination after removal of A2 (which is then orthogonal to every 
N [•] shape). The resulting six shapes are plotted in figure 10. The shapes Ai^ ^3, Aq and 
the mirrored A4, have an overall similarity, while A^ is distinct. 
The mode expansion coefficients, 

of the shape functions, 

Si = N[A^/Kiii, (4.8) 

corresponding to the normalized and divergent- free basic terms in eq. (4.4), are plotted in 
figure 11. The convergence of the expansion series is good (as expected), because we removed 
the divergent local contribution. As the plot shows, O (20) modes are sufficient to achieve 
better than 95% convergence, just a few more as in the case of the equilateral or DBI shape 
(cf. figure 6). Thus, all shapes can be recovered rapidly by the polynomial modal expansion, 
an important requirement for future comparison with data. 

However, the strong correlation, visible to the naked eye, between four shapes in figure 10 
indicates that our ad-hoc choice is not well suited for subsequent discrimination. We can 
make this statement quantitative by computing the correlation coefficients between shapes 
defined via 

C., := , , , (4.9) 

Y {^ii^i)y {Si, Si) 

using the inner product on the tetrapyd in section 2.2.2. The correlation coefficients Cij 
constitute the components of the correlation matrix in figure 12. 



(4.7) 



5 10 15 20 25 

Figure 11. Convergence of the expansion series, i.e. X]n=(r ^ri^ri*/ {^i-> ^i)^ after the N [•] operation 
(removal of divergent terms, substraction of a constant mode and normahzation, cf. eq. (4.5)) on the 
remaining five basic shapes in eq. (4.4). The inset shows that O (20) shapes are sufficient to achieve 
better than 95% convergence. 



Legend 1 ^ ^ 4 5 Legend 




1 2 3 4 5 1 2 3 4 



(a) unordered (b) ordered 

Figure 12. Correlation matrix {Qj}i=i^...^5^ j=i^...^5 of the basic contributions after the N [•] operation 
(removal of divergent terms, subtraction of a constant mode and normalization, cf. eq. (4.5)) on the 
remaining five basic shapes in eq. (4.4). Ordering: (a) ^ 1, ^3 ^ 2, A4, ^3, ... and (b) as 
in eq. (4.10), where ^4 was mirrored, turning correlation into anti-correlation and vice versa for this 
shape. 

Ordering the shapes such that the structure in the correlation matrix becomes more 
visible, namely^^ 

1 ^ N [A] ; 2 ^ TV [A] ; 3 ^ TV [^1] , 4 ^ -TV [^4] , and 5 ^ TV [^3] , (4.10) 

we can identify groups of correlated shapes, figure 12(b), also indicated in the previous line 
by the spacing and the punctuation. 

This grouping agrees with our intuitive eye-balling of the shapes in figure 10. The 
constant mode is, by construction, completely separated. Interestingly, we observe that Aq 
has significantly less correlation with the other shapes, while the correlation among the last 
four terms is at least 78.1%, and among the last three terms above 99.0%. 

As a guide to compare with commonly known shapes, we would like to emphasize that 
the equilateral and orthogonal shapes, 

2 8 

Aqui OC K12 - K3- 2^111 Artho OC Aeqni " ^ Aonst OC Ku - K3 - g^lll (4.11) 

^^One can for example start with an ordering according to the loads of the first eigenvalue and then fine-tune 
the ordering. If extended, this method would lead to a factor analysis or principal component analysis (PCA) 
of the shapes. 
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Figure 13. Plots of the basic normalized and divergent-free contributions N [^4^""^] , where i = 
1, . . . , 6, which are the result of the N [•] operator, cf. eq. (4.5), on the terms in eq. (4.12) and not yet 
plotted in figure 10. These shapes are chosen to minimize the correlation between the basic shapes. 

(see [30], p. 12) are just linear combinations of the N [Ai] and .Aconst shapes. Hence, a better 
choice could involve these two shapes as we already observed their orthogonality. Further, it 
turns out that most of the prefactors of our chosen basic shapes, see eqs. (3.16) and (3.26) 
to (3.31), are slow-roll suppressed. 

Hence, for exactly scale- invariant {ug = 1) single-field models for which the bispectrum is 
built from the basic terms in eq. (4.4), the application of the commonly used shape templates 
provides a good approximation in almost all cases. 

However, being constrained by templates is unsatisfactory for more general/realistic 
models. To elucidate this point, consider the dependent K-terms in eq. (B.4) (without 
considering their potentially slow-roll suppressed prefactors): the correlation with our six 
basic terms is not always high, even though they are a linear combination of the basic 
contributions. This is not surprising, since the basic shapes we chose are not orthogonal. 
Building an orthonormal basis from the latter, as we have done for the separable mode 
functions Qn(x, y, z) (see section 2.2), is clearly an option [24], but the resulting basic shapes 
would be non-separable and increasingly intricate (i.e. after application of Gram-Schmidt to 
the set in eq. (4.10)). 

4.2.2 Finding the least correlated, independent K-terms 

Without invoking an orthogonalization procedure, that is without combining the allowed 
shapes in eq. (4.1), we would like to identify the least correlated, independent K-terms in 
order to estimate their contribution to a given expansion series. We can achieve this goal by 
applying a factor analysis to the full set of possible i^-terms in eq. (4.1); the resulting rotated 
basis 

— Ki2^ A2 — ^const 

.(ic) ^ Ku .(ic) ^ Kg 

and Ks (local shape) is a much better candidate to discriminate observationally between 
shapes without introducing increasingly cumbersome (but fully orthogonal) templates. The 
new terms are A''^''\ A^l''\ and Aq''\ which are plotted in figure 13. Since the new shapes 
remain simple, we find the same rapid convergence for a polynomial mode expansion, see 
figure 14 (only O (20) modes are needed for better than 95% convergence). 

^^Such models are for example the scale- invariant models in [72] and the dominant contributions in slow-roll 
type models as in [16]. 



11, 



(ic) ^ -^123 
^3 ' 



.(Ic) _ — i 

■^3 — 



(4.12) 
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5 10 15 20 25 

Figure 14. Convergence of the expansion series, i.e. X]n=(r ^ri^ri*/ {^i-> ^i)-> after the N [•] operation 
(removal of divergent terms, subtraction of a constant mode and normahzation, cf. eq. (4.5)) on the 
remaining five basic shapes in eq. (4.12), who are chosen to minimize the correlation between the basic 
shapes. The inset shows that O (20) shapes are sufficient to achieve better than 95% convergence. 




1 2 3 4 5 



Figure 15. Correlation matrix {C-^''^}i=i,...,5, j=i,...,5 of the basic contributions after the N [■] oper- 
ation (removal of divergent terms, subtraction of a constant mode and normalization, cf. eq. (4.5)) 
on the remaining five basic shapes in eq. (4.12), ordered as in eq. (4.14). These shapes are chosen to 
minimize the correlation between the basic shapes. 

The correlation matrix 

gg':= , ^t^Sp ^here 5^ = iV [^H /^ni' (4-13) 

for the ordering 

1 o -N [Af^] [A'i''>] , 4 o iV [4'='] , and 5 o -A^ [^^'^'j , 

(4.14) 

is plotted in figure 15: correlation between 62.4% and 88.1% on the off-diagonal and 81.0% 
anti-correlation between N [^j^] and N [i^i2] is evident, a significant improvement compared 
to figure 12. 

We would like to highlight the usefulness of such a least correlated, but simple choice 
of basic shapes. A general (scale-invariant) bispectrum is a linear superposition of the basic 
shapes 

A^'^Ci N [Af"^] + Cconst Aonst + Qocal v4local- (4.15) 

i 

By construction, model parameters, such as 05,65,..., enter only via the coefficients q, 
Cconst and Qocal 5 ^ot in Af''''. Let's assume we would like to identify the most likely model 
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parameters given a measurement of modal expansion parameters. To scan the parameter 
space, e.g. via MCMC simulations, modal expansions for many As are needed; fortunately, 
this requires a one-time modal expansion of the basic shapes only. At this step, any choice 
of simple basic shapes suffices, but if the chosen basic shapes are only weakly correlated, 
a resulting bispectrum is less likely the result of a near perfect cancellation of some highly 
correlated shapes (the latter would amplify errors). Furthermore, it is possible to pinpoint 
the basic shape that dominates the bispectrum, enabling physical insight into the origin 
of non-Gaussianities; the latter would be obfuscated if a cumbersome orthogonal set were 
chosen. 

Roughly speaking, the least correlated basic shapes provide the middle ground between 
the orthogonal mode functions on the tetrapyd, whose expansion coefficients are optimally 
constrained by the data without, however, offering direct physical insight, and the full primor- 
dial shape function, which is the direct imprint of a cosmological model, but unfortunately 
not directly accessible in observations. 

4.3 Expanding slow-roll violating models 

Equipped with a thorough understanding of general shapes in single-field slow-roll models, 
we come to the main focus of this paper, the slow-roll violating single-field models introduced 
in section 3.3. We are interested in models with a bispectrum that is a mixture of different 
commonly discussed shapes (cf. [23], figure 2). To this end, the discussion of shapes in 
section 4.2.2 is directly applicable as we use the least correlated basic shapes for an optimal 
reconstruction of the shape coefficients using the modal decomposition. We constrain our 
investigation to models that are consistent with the newest WMAP data, that is we assume 
an almost scale-invariant power spectrum with a slight red tilt. The WMAP7 result for the 
spectral index is [19] 

Us ^ 0.968 lb 0.012 (68% CL) (WMAP7 only). (4.16) 

As concrete examples, we would like to focus on two models: the general DBI model 
with f^^'-' = 1 — c^, and the negative fx models where we chose fx = —100 as in [23], 
figure 4. In both models, we chose a small speed of sound, namely Cg = 1/20, in order to 
have a non-negligible f^^. The slow-roll parameter e is chosen in the range [0,0.3], so that 
observations remain in line with CMBR observations (cf. [23], eq. 2.9). 



Model 


€ 


Cs 


fx 


Us 


General DBI model 


{0.3, 0.2, 0.1, 0.075, 0.05, 0.04, 0.03, 
0.025, 0.02, 0.015, 0.01, 0.005, 0.001} 


0.05 




0.96 


Negative fx model 


{0.3,0.15,0.060,0.03, 
0.015,0.006,0.003,0.001} 


0.05 


-100 


0.96 



Table 2. Parameter table for two exemplary slow-roll violating models to which we apply the modal 
expansion. 

Our first model corresponds to the first concrete example chosen in [23], section 4.1 
(fast-roll DBI inflation). For models in this class, the first term in the bispectrum, ^^3, is 
negligible compared to the other amplitudes in eq. (3.39). 

The negative fx model is an example from the class of models discussed in [23], sec- 
tion 4.2 (|A/5]| ^ 1). The equilateral shape is the most prominent shape in this class. 
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equi 

.equi Jnl 




(a) general DBI model (b) negative fx model 

Figure 16. for the two concrete examples of the Noller and Magueijo shape, (a) The general 

DBI model with = 0.96, Cs = 1/20 and /j^"^'^ = 1-C5^ (b) The negative fx model with = 0.96, 
Cs = 1/20 and fx = -100. 




Figure 17. Plots of the basic normalized and divergent-free contributions L [Af""^] , where z = 1, . . . , 6, 
which are the result of the L [■] operator, defined beneath eq. (4.17), on the terms in eq. (4.12) except 
-4iocai = These shapes are chosen to minimize the correlation between the six basic shapes. (Note 
the different scales of the vertical axis.) 

dominating over the remaining ones such as the local contribution caused by the violation of 
scale invariance. 

In both cases deviations from slow roll decrease the magnitude of see figure 16.^^ 

4.3.1 Approximative scale invariance 

In the limit of exact scale invariance, the Noller and Magueijo bispectrum reduces to the one 
by Khoury and Piazza [72]. In this limit, the {Kiii/2K^)'^^~^ prefactor vanishes for all shape 
functions in eq. (3.39). Hence, the full bispectrum ^nm is a linear superposition of the base 
K terms with coefficients depending on the parameter set {e, Cg, fx}- 

^^We employ, just as in section 3 and in [23], the /nl sign convention by Komatsu et al. [19]. 
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Figure 18. Convergence of the expansion series, i.e. X^n^^ ^n^^n^*/(^f°^^ ^i^"^^)? terms in 

eq. (4.17). The inset shows that O (20) shapes are sufficient to achieve better than 95% convergence. 



Since we are discussing the observationally motivated case of a small deviation from 
scale invariance, one can expand the full bispectrum around cti = - 1 = 0.2° The leading 
order contribution is the scale-invariant case discussed in [72], whose basic terms are given in 
eq. (4.4). To subleading order, i.e. O (o^i), a linear combination of the basic K terms and an 
additional set of base shapes that are the product of the ones in eq. (4.4) and log(Kiii/K^) 
arise. Such terms are also present in the bispectrum by Chen et al. [16]. 

Because the multiplication with the logarithm should not drastically affect the overall 
shape, we expect that these new terms are strongly correlated with their "mother terms". 
In order to quantify this conjecture, we employ the set of least correlated, independent basic 
terms in eq. (4.12) and calculate the correlation matrix for these terms and the corresponding 
additional contributions. 



L [K12] , L [Km] , L 



K 



L 



Kl4 



L 



Kg 



and L 



Ku3 



(4.17) 



here the operator L [•] removes the local contribution from a given shape, multiplies the 
remainder with log(A^iii/A^^) and normalizes the product using the N [•] operator. Note 
that f^i^^^ is non-divergent for all terms in eq. (4.17). The convergence of the expansion 
series of the shape functions 

Sf°^ = L [^f^^] /i^ni, (4.18) 

is again rapid, figure 18, and we achieve better than 95% convergence using 20 mode functions 
for all terms in eq. (4.17). 

The full bispectrum up to first order in oli contains only the terms in Eqs. (4.12) and 
(4.17), i.e. 

-4,nm(01, 0:2, Cs, fx) = /local Aocal + /const ^const 

+ E./.iV[^n+E.5.^K1+0(a?), 



(4.19) 



where the /... and gi are of order O (ai). The S^div — K^log^Km/ K"^) shape, whose f^f^^^ is 
not defined because \og{Kiii/ K^) cx), is absent. ^-'^ Hence, /^l^^ of the whole expression 

^"^Note that the expansion coefficients o^n, see eq. (2.19), (and all derived variations such as an and ol^^) 
should not be confused with the model parameters oli and 0L2 in eq. (3.36). Aiming at full consistency with 
the literature, we cannot avoid this conflicting notation. We distinguish the model parameters ai and CX2 by 
a bold face in order to guard against misunderstandings. 

^^Naively, there is a formal \og{Kiii / K^) term in the O (ai) bispectrum which is, however, canceled by 
the L [•] operator. 



1 2 4 6 ^ \° Legend 




12 4 6 8 10 11 



Figure 19. Correlation matrix {C-^°^^}i=i,...,io, j=i,...,io of the terms in eq. (4.17), ordered as in 
eq. (4.20). 



up to first order in ai is proportional to /locah which can be derived from the consistency 
relation in eq. (3.17) by Creminelli and Zaldarriaga [74] (i.e. the local contribution is absent 
in the case of exact scale invariance). In the following, it is important to keep in mind that 
even though the almost scale-invariant case leads to a local contribution, it is suppressed by 
the smallness of — 1. 

The correlation matrix in figure 19 is best visualized if we order shapes as 

1 ^ -N [Ku] ; 2 ^ -TV [Ku/K^] , 3 ^ -L [Ku/K^] ; 
4 ^ L [Km] , 5 ^ L [Kus/K^] , 6 ^ TV [Kus/K^] ; .Q^ 
7^L [Kq/K^] , 8 ^ L [K^/K] , 9 ^ TV [Kq/K^] ; ^ ' 

10 ^ L [Ku] , 11 ^ TV [K4/K] , 

where correlated groups are indicated by spacing and punctuation. The structure of the 
correlation matrix proves, as conjectured, the similarity of the additional shapes to their 
mother terms (or another basic term). 

For example, TV [Ku/K'^] and L [Ki^/K'^] have a 99.9% correlation, but there are three 
exceptions from the general rule: the shapes L [A^m], L [K12] and L [K4/K]. 

^[^111] is, by definition, orthogonal to Km. This shape is 99.8% correlated with 
L [Ki23/K^] , which itself is only 96.1% correlated with its mother term TV [^123/^^^] . This 
relationship is apparent in the plots (compare figure 17(b, f) and figure 13(c)). 

L[K4/K] shares a higher correlation (96.2%) with TV [Kq/K^] than with L [Kq/K^] 
(89.8%). This is related to the behavior of the shapes in the squeezed limit (compare fig- 
ure 17(c, e) and figure 13(a)). 

The mother term of the former, N[K4/K], is 96.9% correlated with L[Ki2], because 
the divergence of S in the squeezed limit is suppressed. The latter is only 77.4% correlated 
with its mother term, TV [i^i2]- 

4.3.2 Convergence of the NoUer and Magueijo example models 

Because observations show a nearly scale-invariant power spectrum eq. (4.16), we restrict 
our analysis to models that satisfy this constraint, taking Ug = 0.96 as a reference value. 
Consequently, the {Km/ K^)'^^~^ factor in eq. (3.39) leads to a divergent contribution to 
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(a) General DBl expansion 
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(b) General DBl convergence 
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(c) Negative fx expansion 



(d) Negative fx convergence 



Figure 20. Coefficients and convergence of the expansion series, i.e. X]^=o'' ^n^^n^*/(^i5 ^i)? 
bispectra for the models defined in section 4.3 for selected model parameters. 



the bispectrum, which cannot be subtracted easily and we have to fall back to the clipping 
method introduced in section 2.3.^^ 

We calculated the expansion series of the two proposed example models, the general 
DBl model and the so-called negative fx model, for the set of parameters given in table 2. 
The expansion coefficients of the normalized shape functions, 

5'NM(fci, ^2, ^3) = where S^uik, k, k) = 1, (4.21) 

N [kik2k^Y 

in the case of the general DBl model are given in figure 20(a) for selected values of the slow 
roll parameter e; they are compared to the usual DBl model and the clipped local model 
(divided by a factor of 3 in order to have compatible amplitudes). 

In the same manner, we give the expansion coefficients of the negative fx model in 
figure 20(c), compared to the equilateral and the clipped local model (the latter is again 
divided by a factor of 3). 

In line with the dependence of the equilateral non-linearity parameter f^^ on e (see 
figure 16), we find that the relative equilateral contribution decreases with increasing e. The 
local contribution is almost constant because the spectral index is unchanged. We did not 
include error bars in the plots of the expansion coefficients in figure 20(a, c), which are given 
by the accuracy of Mathematica's numerical integration utility. In order to have physically 

^^Substraction was possible in the last section, since the bispectrum had been expanded to first order in ai 
only; here, we would like to perform the modal expansion of the full bispectrum, that is, to go beyond first 
order in cki. 
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(a) base term contribution 



Si 




(c) reconstructed contributions 



Figure 21. (a, b) Contributions of the terms in eq. (4.19) to the bispectrum amphtudes for the 
general DBI model defined in section 4.3 for a certain set of example parameters, (c) Results of the 
linear regression on the expansion coefficients as in eq. (4.23). 

significant error bars for the mode coefficients relevant for CMB data comparison, one could 
derive the uncertainty from simulated maps [27]. 

The relative increase of the local contribution explains why the correlation of the ex- 
pansion series with the full bispectrum, shown in figure 20(b, d), decreases for larger e. For 
the largest deviation from slow roll, we still achieve better than 92% convergence using 32 
modes for both example models; we achieve almost 93% using the full rimax = 53 modes (not 
shown) . 
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Figure 22. Calculated, predicted (up to O (cxi)) and recovered (from a linear regression using 
reference shapes) expansion coefficients of the general DBI model with e = 0.075. We plot the 
coefficients of the bispectrum amplitude (not normalized) in order to compare with predictions in the 
range n = 13, . . . , 33. 



4.3.3 Consistency of the relative contributions of the basic terms to the full 
bispectrum 

In order to compare the expansion coefficients with our expectations, we expand the full 
bispectrum as given in eq. (3.39) around ai = — 1. We decompose the bispectrum 
amplitudes of order O (ai) into the least correlated basic terms and the derived L [-J-type 
shapes as in eq. (4.19) in order to improve the reconstruction that we employ to verify our 
results. 

In figure 21 we plot the contributions of the basic terms to the zeroth and first order 
bispectra for the discussed set of model parameters. 

During the discussion of the L [-J-type shapes in section 4.3.1, we learned that we cannot 
distinguish certain shapes due to their high correlation. We thus take one representative 
from each group and estimate the contribution of representatives by a linear regression on 
the expansion coefficients. The result of that procedure can be seen in figure 21(c). The 
dashed line is the expected combined contribution for each group and the full lines are the 
results of the regression on the fi parameters, 

arg min (^an - /i c^^I^ j Vn < n^^ax = 52, (4.22) 

where the are the expansion coefficients of the shapes representing a group of highly 
correlated shapes (compare figure 19); representatives of our "binned" shape classes are: 
the constant mode, the clipped local template, N N [K^]^ N [Ku]^ N [Kq/K^]^ and 

N [K222/K^]- The numerical inversion of the correlation matrix would be unpredictable if 
we had chosen to reconstruct the individual contributions of the highly correlated shapes. 
The result of this regression yields^^ 

'^shapes 

/i= E where (4.23) 

i=i 

^max 

= ^n^^n Vz, j = 1, . . . , rishapes and (4.24) 

n=l 

^^We employ no weights for the coefficients since we assume them to have comparable errors. 
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fi 




Figure 23. (a,b) Contributions of the terms in eq. (4.19) to the bispectrum amphtudes for the 
negative fx model defined in section 4.3 for a certain set of example parameters, (c) Results of the 
linear regression on the expansion coefficients an as in eq. (4.23). 



■^max 

dj = ^ (^n^Oln Vj = 1, . . . , nshapes- (4.25) 

n=l 

The reconstructed coefficients deviate slightly from the predictions as shown in fig- 
ure 21(c), but the order of the contributions has been fully recovered. The two largest 
contributions, the constant mode and the N [K12] shape (corresponding to the equilateral 
and orthogonal template) are well reconstructed. The 95% confidence intervals (assuming 
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Figure 24. Calculated, predicted (up to O (cxi)) and recovered (from a linear regression using 
reference shapes) expansion coefficients of the negative fx model with e = 0.006. We plot the 
coefficients of the bispectrum amplitude (not normalized) in order to compare with predictions in the 
range n = 13, . . . , 33. 

Gaussian distributed an) of the linear regression results are indicated by the shaded regions. 
Corresponding to the large correlations among the base shapes, the confidence intervals are of 
the same order as the reconstructed values for all shapes except for the (orthogonal) constant 
mode, the (clipped) local template, and the N[Ki2] shape, probably caused by a systematic 
effect due to the employed clipping mechanism. 

Why are the N [Ki4/K'^] and N [Kq/K^] contributions not well recovered? There are 
systematic deviations for both contributions that can be explained as follows: first of all, 
our predictions are first order in oti only, which can lead to small deviations. Second, we 
neglect the differences between the "binned" and referenced shapes. We find that this effect 
is insignificant for the general DBI model, as the expected mode coefficients and the sum over 
binned amplitudes as well as reference mode coefficients are almost the same, see figure 22 
for an exemplary reconstruction for e = 0.075. However, the differences between coefficients 
of the O (cxi) prediction and the original expansion coefficients can explain the systematic 
deviation of the reconstructed contributions. 

We likewise reconstructed the bispectrum amplitudes for the negative fx model. Again, 
the order of the contributions is fully recovered and the two largest contributions are well 
recovered quantitatively (see figure 23). The local contribution is, as expected, relatively 
small compared to the equilateral one. 

We can explain the large confidence intervals for the base shapes with the same reasoning 
as for the general DBI mode, that is by large correlations. Systematic deviations are present 
for the N [Ki^/K'^] and N [Kq/K^] contributions, caused again by the same systematic 
effects explained above. Looking at the calculated, expected and reconstructed coefficients 
for the exemplary case e = 0.006 (see figure 24), we see that binning does not alter results 
significantly, but the expansion coefficients of the full bispectrum amplitude deviate from the 
O (oti) expectation, leading again to a systematic bias in the reconstruction. 

We conclude that the reconstruction of dominant contributions in expanded bispectra, 
the standard templates in the case at hand, is not only qualitatively, but also quantitatively 
possible. The reconstruction of the remaining subleading basic shapes is more challenging. 
The main drawback we encountered is the large correlation between basic shapes, which 

^^We use Mathematica's LinearFitModel. The confidence intervals of the parameters are computed using 
student-t statistics, where the parameter's standard errors are taken from the diagonal terms of the covariance 
matrix. 
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weakens constraints; furthermore, we found systematics in the not exactly scale-invariant 
cases, which should be incorporated into the confidence intervals. 

4.4 Summary of our results 

We first identified the basic contributions to the bispectrum amplitude A of scale-invariant 
fast-roll models, which are of the form i^3, Kpq/KP^^~^ and K^st/ • These shapes are 
degenerate (linearly dependent) so that their relative amplitude in an observed bispectrum 
cannot be separated. We found that naively chosen linearly independent basic terms are 
prone to high correlations: in figure 12(b), it is clear that three terms in eq. (4.4) are of the 
equilateral type and can hardly be distinguished from the equilateral template. In this case, 
we are essentially left with the three common templates (local, equilateral and orthogonal) 
that measure the divergent contribution removed from our N [i^...] shapes, the equilateral- 
type shapes and the constant mode. Thus, the contribution of Aq — ^222/^^ is the only one 
left to be analyzed after measuring the three common /nls, because A^ does not contribute 
to the bispectrum amplitude in eq. (3.39). 

Fortunately, this predicament can be alleviated considerably: to optimize observability 
and separability of the shape constituents, we identified the least correlated basic shapes 
(factor analysis) that spans the same space of allowed shapes, see eq. (B.4). This step avoids 
block structures in the correlation matrix (see figure 15), leading to an invertible matrix. 
As a consequence, we can estimate the amplitudes of the local template, the constant mode 
and all five remaining independent shape constituents in a set of expansion coefficients by 
a linear regression. Furthermore, as these least correlated shapes are also well described by 
a truncated expansion series (rapid convergence), one can reconstruct expansion coefficients 
of general bispectra in scale- invariant single-field models fast, which is crucial for efficient 
MCMC simulations 

Shifting our focus to almost scale- invariant models, we took into account first order 
corrections in n^ — 1, leading to subleading contributions in the bispectrum of the L [i^...] form. 
To first order, we show that these subleading contributions fall into subclasses composed of 
the least correlated basic shape constituents, as these terms are strongly correlated with 
either their "mother terms" or another basic shape. 

The convergence of the expansion for our example models from table 2 is challenged 
by /nl^H-4nm) = 00 (unless only terms up to first order in O [oci) are taken into account), 
which required the use of a clipping technique. The correlation of the expansion series and 
the full bispectrum plotted in figures 20(b) and 20(d) shows that the more relevant the local 
contribution is, the lower the degree of convergence becomes. 

Nevertheless, we showed that reasonably good convergence can still be achieved for two 
exemplary model families (not exactly scale-invariant, not fine tuned for convergence in any 
way): they can be described by (9(50) expansion coefficients with a loss of information of 
at most 10% - even the largest deviation from slow roll (having the highest relative local 
contribution) is 90% correlated with the full bispectrum while truncating the expansion 
series at rimax = 53. We find that contributions of the dominant constant mode, the N [K12] 
shape, and (to a lesser degree) the local template can be recovered quantitatively, while a 
reconstruction is less successful for the subdominant basic shapes (even if the least correlated 
basic shapes are used). 

The dominant contributions are exactly the ones that are well described by the projected 
non-linearity parameters /^l^^' /nl^^5 hence, for this particular class of single- 

field models it is sufficient to focus on these three well known non-linearity parameters. The 
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correlation between subclasses of the remaining shapes leads to large uncertainties on the 
best-fit parameter. In addition, we find systematic effects complicating reconstruction. 

To summarize, we provided a complete dissection of possible bispectra in nearly scale- 
invariant, single-field models of inflation. An important tool are modal techniques, which 
were used systematically. Further, we would like to stress the importance of using the least 
correlated basic shape constituents as an intermediate step, both for subsequent efficient 
MCMC simulations as well as for reconstruction. The toolchain employed in this section is 
completely general and can be applied to more complicated bispectra (primordial, CMB or 
late time). 

4.5 Relationship of our analysis to recent work 

While completing this paper, a related article [24] appeared. Ribeiro and Seery discuss the 
phenomenology of non-Gaussianities of general scalar field models, whose curved background 
extensions maintain second-order field equations and stress tensors. The corresponding La- 
grangian was first written down by Horndeski [69] and recently resurfaced within the so-called 
Galileon inflationary models where a Galilean shift symmetry is required. This model gained 
popularity over the last years [70] (see [54] for related work). The possible NG bispectrum 
shapes for general Horndeski models are identified in [80-82]. 

Ribeiro and Seery [24] relate a newly found orthogonal shape to previously described 
orthogonal models — orthogonal to the equilateral and the so-called orthogonal template — by 
means of a mode decomposition; hence no qualitatively new non-Gaussianity shapes appear in 
these models. The cases discussed in [24] are related but not identical to the ones investigated 
here. Furthermore, our main focus has been the systematic dissection of bispectra. Ribeiro 
and Seery emphasize a particular new shape that appears in their models and investigate 
similarities to related shapes. 

4.6 Outlook to application in N-body simulations 

Non-Gaussianities are present in large scale structure surveys, but it is challenging to disen- 
tangle late-time non-linearities from primordial signals. Nevertheless, observational bounds 
based on the scale dependent bias [5-8] are becoming competitive with CMB results; to 
improve upon these bounds, N-body simulations are needed to shed light onto the dynam- 
ical evolution of non-Gaussianities. Current simulations use primarily non-Gaussianities of 
the local type to set the initial conditions for the gravitational potential. However, if non- 
Gaussianities are observed, it becomes of prime interest to discriminate between different 
primordial shapes and thus be able to invoke more general initial conditions. 

Just as for CMB physics, factorizable templates (local, equilateral and orthogonal) lead 
to reduced computational effort. However, for more general shapes brute force methods 
appeared to be the only option [83, 84] for a while. Given the tools employed in this paper, 
it is not surprising that one can remove this hurdle by applying modal techniques, as shown 
recently in [39].^^ 

We are currently implementing [40] modal techniques for the bispectrum as well as 
the trispectrum and aim to run large simulations in the foreseeable future with the goal to 
follow the evolution of shapes as non-linearities grow, ultimately leading to the use of N-body 
simulations in conjunction with LSS surveys to constrain primordial non-Gaussianities and 
thus inflationary models. 

^^The kernel used in [39], which affects the modal decomposition, is only one possibility; further details are 
given in [38]. 
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(a) real part 



(b) imaginary part 



Figure 25. Plot of the first ten orthonormal one-dimensional Fourier modes, Qr (see [37], figure 6) 
using a Fourier standard basis instead of monomials. 

5 Oscillatory and other features in the bispectrum, alternative mode func- 



In Section 2.4, a specific scale-dependent model has been introduced, the sharp feature model. 
Due to slow convergence w.r.t. the polynomial mode functions, see figure 9(b), a set of 
alternative mode functions (Fourier modes) was proposed by Meerburg [37] . 

In this section, we review Meerburg's choice of Fourier modes and explain why a selected 
subset of the mode functions increases the correlation between the truncated expansion series 
and the full amplitude. Furthermore, we apply the alternative modal expansion to a recently 
discovered shape that stems from particle production during inflation. 

5.1 Fourier mode functions 

After reviewing the construction of the alternative mode functions, we discuss the improved 
convergence for oscillatory bispectrum features. 

5.1.1 Construction of the orthonormal basis 

Using Fourier modes as in [37], the construction of the orthonormal basis in section 2.1 needs 
some minor modifications. First, we replace the monomial standard basis in eq. (2.7) by a 
Fourier standard basis. 



in the definition of the orthonormal one-dimensional functions (orthonormal w.r.t. the weight 
function w). Since the mode functions become complex, the second argument in all inner 
products (/, g) is complex-conjugated. 

Denoting all functions based on these Fourier modes by a tilde, we plot the tetrahedal 
Fourier modes qr{x) for r = 1^ ... ^10 in figure 25. The next steps are in complete analogy 
to the polynomial case in section 2.1. We apply the symmetrization procedure to compute 
the Qprs{x,y, z) and Gram-Schmidt orthonormalization in order to obtain the orthonormal 
Fourier mode functions Rn{x,y,z) on the tetrapyd domain. 

5.1.2 Convergence properties of common shapes 

Because we expand a real shape function S in complex modes, the modal coefficients 



tions 




r G No 



(5.1) 



{Rn, S) 
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Re(a„), Im (an) 
0.10 r 



equilateral shape [real part] 
equilateral shape [imaginary part 
DBl model [real part] 
DBl model [imaginary part] 




Figure 26. Plot of the expansion coefficients for the equilateral shape and the DBl model w.r.t. the 
Fourier mode functions Rn. Note that the constant mode contributions lie outside the plot range at 
ai = 0.424 and ai = 0.373, respectively. 




equilateral shape 
DBl model 



(a) equilateral shape and DBl model 




Figure 27. Plot of the sum of squared coefficients (a) for the equilateral shape, see eq. (2.21), and 
the DBl model, see eq. (2.22), and (b) for the feature model as given in eq. (2.34) for = k^a^/ujf^ 
where cj/ G {8, 9, 10, 15, 20, 30}, and 6 = 0. Both expansions are w.r.t. to the Fourier mode functions. 



fulfill 5* = a-n-'^^ Parseval's theorem (see eq. (2.24) for the polynomial case) needs a slight 
modification, i.e. 

oo 

(5, 5) = 5? + 2 ^ Re [5^(5;] . (5.3) 

n=2 

Due to the oscillatory nature of the Fourier modes, we cannot expect a good conver- 
gence of our highly correlated test shapes, the equilateral / DEI shapes in eqs. (2.21) and 

^^This equality is the reason, why we still can take r > 1 and n > 1. 
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(2.22). In fact, apart from the constant mode coefficient ai (accounting for 81.7% and 89.4%, 
respectively), the first 53 Fourier modes do not contribute much, as shown in figure 26. We 
find 91.6% and 96.1% correlation, respectively, between the truncated expansions (nmax = 53 
modes) and the respective full shape. 

Hence, the degree of convergence is worse compared to the polynomial expansion for 
these shapes, as shown in figure 27(a). In terms of the convergence indicator s introduced 
in section 2.2.1, we find s^^^i^^^ ^ 0.084 and ^^^^^^ ^ 0.040 after rimax = 53 modes, clearly less 
satisfactory than the polynomial expansion. 

However, the replacement of polynomials by Fourier modes pays off in the case of 
localized and/or oscillatory contributions to the bispectrum, such as in feature models [32, 85] 
and resonant non-Gaussianities [32, 33, 64, 65]. We agree with Meerburg [37] regarding the 
improved convergence properties, e.g., for the sharp feature models defined in eq. (2.34). We 
plot the convergence of the Fourier mode expansion series for the cases fc* = fcmax/^/, where 
(jjf G {10, 20, 30, 40}, in figure 27(b). In comparison to the polynomial expansion in figure 9, 
we see that particular modes improve convergence drastically, while the majority do not. As 
Meerburg pointed out, this is due to the fact that the sharp feature model is analytically 
expandable in separable trigonometric functions by using trigonometric identities (eq. 22 in 
[37]). For lower frequencies, the convergence is improved compared to the polynomial case 
shown in figure 9(b). For higher frequencies both expansions are futile, since the oscillations 
are too rapid to be recovered by the first mode functions. 

Meerburg investigated two other classes of models. The first group encompasses the 
resonant-type non-Gaussianities that arise due to a periodic feature in the inflaton potential 
[32, 33, 64, 65] (present, e.g., in monodromy inflation [64, 86]). He took 



as an example shape and found good convergence for the lower part of the allowed frequency 
range (20 < cj^ < 10^). The number of modes needed is decreased by a factor of 5. However, 
he did not achieve more than 15% convergence using 80 Fourier modes for uj^ = 80, with 
worse results for higher frequencies. As a consequence, non-Gaussian signals of the resonant 
type might not be detected within an analysis that is constrained to (9(100) modes, even if 
Fourier modes are used. 

The second group are non-Gaussian shapes that arise in theories with initial state mod- 
ifications [16, 53, 58, 66] where small deviations from the Bunch Davies vacuum can lead to 
large NG signals. The reconstruction of such models is challenging with either basis, because 
a large number of features coupled with rapid oscillations appear. Furthermore, features 
in different regions of the tetrapyd have different effective frequencies, introducing a scale 
dependence (see [37], figure 1). 

5.1.3 Pre-ordering mode functions 

For resonant bispectra of the type in eq. (5.4) (not analytically separable via trigonometric 
identities) only a shape-specific fraction of mode numbers (see figure 7 in [37]) is relevant. For 
this particular shape, the relevant modes are those that oscillate along the K = fci + A:2 + ^3 
direction, consistent with the analytic expression, see eq. (5.4). 

Thus, if a shape is restrained to some specific form, one may simply pre-order the 
mode functions according to their relative contribution to the overall convergence so that 
considerable less modes are needed for the same degree of convergence. 




(5.4) 
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Maybe more importantly is the application of such a technique to the actual data, 
such as the CMB bispectrum, which contains known acoustic oscillations. Applying a factor 
analysis to the PLANCK data in order to identify the important modes in a polynomial 
expansion is currently being considered. 

5.2 Particle-production scenarios 

We would like to extend the analysis of oscillating, scale-dependent shapes to signals that arise 
due to particle-production scenarios during inflation. Because of interactions with additional 
fields, non-inflaton particles can be produced at special points on the inflaton trajectory. The 
backscattering of these particles leaves a trail in the bispectrum of curvature perturbations 
in the form of localized features (see [26] for a review). Such events are common in trapped 
inflation [43, 44, 46, 87], monodromy inflation [48], or modulated perturbations from extra 
species point encounters [45, 47, 88]. 

In order to examine those scale-dependent bispectra, we focus on a shape derived by 
Barnaby et al. [41, 42] for the case of a simple prototype interaction 

Ant = -y(</'-</'o)V (5.5) 

between the inflaton, 0, and an iso-inflaton, x (see [89] for a more realistic model, containing 
interactions of a pseudo-scalar inflaton with gauge fields). The shapes of the resulting non- 
Gaussianities are discussed in [25] (see [90] for a recent moment analysis). The imprints in 
the primordial bispectrum share qualitative features with the ansatz 



^ 1 — cos A / + rn^t 

-7r/ef/(3fc2)_ 



kf + m? 



(5.6) 



where k^ defines the characteristic scale at which the feature occurs. The interaction takes 
place near the point (f) = (f)o, so that we can take 

0(t) ^ 00 + (5.7) 

Following [25], the characteristic scale is defined as k^ = ^/oWl^ so that fc* ^ 30H in the 
case ^ 0.1. Assuming a standard chaotic inflation model, ¥{(/)) = m^(/)^/2 with m ^ 
10~^y/87TMp and (f) = 3.2\/87rMp, Barnaby used lattice simulations to extract the probability 
density function (PDF) for Sep from the simulated data. In addition, he was able to calculate 
the bispectrum analytically. Eq. (5.6) is a simplified expression that resembles the qualitative 
behavior, but should not be used for comparison with data. Backscattering also affects the 
power spectrum, leading to a "bump feature" at a wavenumber fcbump specific to the location 
of the interaction on the inflaton trajectory. Depending on the relation of the wavenumber 
to the observable window, the shape of the contribution to the bispectrum differs (figure 5 
in [25]). The feature in the shape function is localized in any case, and oscillatory for 
fcbump < kma^ (lower right plot of figure 5 in [25]), corresponding to fc^ax ^ e^fcbump)- 

Taking 3.2V87vMp and m IQ-^VSnMp, while using g'^ 0.01 and A:* 30i7, 

we analyzed the shape in eq. (5.6) at the time t — jj and at the location fcmax = e^fcbump 
with our mode expansion tool-chain and tested it for convergence, both with tetrahedral 
polynomials and oscillating Fourier modes. 



^The full, rather cumbersome bispectrum is also derived in [25], eq. 81. 
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1.0^ 




(a) shape (b) convergence 



Figure 28. Plot of the sum of squared coefficients for the particle production model as proposed in 
[25]. 

It is evident from figure 28 that both expansions achieve a higher than 85% convergence 
using rimax = 53 modes. This could be expected, because the ansatz in eq. (5.6) is already 
separable and of oscillatory form. This explains the sharp rise in figure 28(b) at the fourth 
Fourier mode; the exponential damping along the K = fci + /c2 + ^3 direction prevents a 
faster convergence of the expansion. The polynomial expansion converges slower but reaches 
almost the same level using rimax = 53. 

We would like to emphasize that the Fourier expansion is able to tell apart the sharp 
feature model from the particle production scenario by the indices of the relevant modes: 
these are the modes along the K direction, like the fifth mode, for the feature models, 
whereas our specific choice of a particle production scenario peaks at the fourth mode (which 
is (0, 0, 2), see appendix A). As argued in section 5.1.3, a good expansion should be re-ordered, 
assigning lower indices to modes that are most likely relevant for a particular shape family. 
These relevant mode coefficients should then be determined from the data and compared to 
the predictions. 

We plan to perform a more detailed investigation of the full bispectrum resulting from 
both singular or repeated particle production events during inflation in a future publication. 
In that regard, the relatively good convergence of the Fourier expansions is encouraging. 

6 Conclusion 

Confronted with increasingly involved models of the early universe and equally complex non- 
Gaussianities, i.e. the bispectrum shape, it is a daunting task to tell models apart based on 
limited observations of the cosmic microwave background radiation or large scale structure. 

We provided an in-depth investigation of a complex, non-separable primordial shape 
family (in fc-space) originating in general single-field models of inflation (Lorentz invariant 
Lagrangian, at most first order in field derivatives [23], reviewed in this paper). We disen- 
tangled the general bispectrum shape into a minimal set of simple basic shape constituents 
(independent of model parameters), chosen to minimize correlation among each other by 
means of a factor analysis. These shapes show rapid convergence in a polynomial modal 
expansion on the tetrapyd domain (expansion techniques of [27] are reviewed in this paper). 
Once expanded, the corresponding modal coefficients of a general shape can be computed 
efficiently via a simple sum, a great advantage for subsequent MCMC simulations. Fur- 
thermore, the basic shape constituents offer physical insights into the general shape, since 
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they can (in principle) be discriminated in observations. Roughly, the least correlated ba- 
sic shapes interpolate between the orthogonal mode functions (optimally constrained by the 
data, but opaque w.r.t. physical processes) and the full primordial shape function (imprint 
of a cosmological model but not directly accessible ). 

We investigated two concrete, slow-roll violating model families (general DBI and fx 
models [23]) resulting in large non-Gaussianities. We provided modal expansions of these 
shapes and identified the dominant basic shape contributions. We find new shapes that 
are, however, closely correlated to known templates or suppressed for the models under 
consideration. Hence, in this concrete class a restriction to common separable templates 
(local, orthogonal, equilateral) provides a good approximation. The techniques demonstrated 
in this case study are readily applied to other, more complex models (in preparation) and 
can also be used for late-time shapes in multipole space, as needed for comparison with the 
CMB data. 

In addition, we briefly summarize current work on N-body simulations, for which we 
employ the same modal techniques. 

However, not all primordial shapes show the needed rapid convergence in a polynomial 
modal expansions. As a concrete, hitherto untreatable case we discuss the bispectrum orig- 
inating from particle production during inflation (more precisely caused by the associated 
backscattering of produced particles from the inflaton condensate), as present in trapped 
inflation or monodromy inflation. This shape contains oscillatory as well as localized fea- 
tures. These can be recovered if Fourier modes are chosen instead of polynomials, which we 
show explicitly. Further, model discrimination is simpler, since key frequencies are readily 
identified in a Fourier expansion. Hence, it might be prudent to perform future data analysis 
with more than one basis set that complement each other, such as polynomials and Fourier 
modes. In addition, proper ordering of modes w.r.t. their relative importance to the data 
could improve convergence and thus efficiency. 
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A Slicing ordering of the mode functions 

We use the so-called slicing ordering of the mode functions Qn = Qprs'i up to nmax=53 it is 
(cf. [27], eq. 58) 

1 ^ {0, 0, 0}, 2 ^ {0, 0, 1}, 3 ^ {0, 1, 1}, 4 ^ {0, 0, 2}, 5 ^ {1, 1, 1}, 

6 ^ {0, 1, 2}, 7 ^ {0, 0, 3}, 8 ^ {1, 1, 2}, 9 ^ {0, 2, 2}, 10 ^ {0, 1, 3}, 

11 ^ {0, 0, 4}, 12 ^ {1, 2, 2}, 13 ^ {1, 1, 3}, 14 ^ {0, 2, 3}, 15 ^ {0, 1, 4}, 

16 ^ {0, 0, 5}, 17 ^ {2, 2, 2}, 18 ^ {1, 2, 3}, 19 ^ {0, 3, 3}, 20 ^ {1, 1, 4}, 

21 ^ {0, 2, 4}, 22 ^ {0, 1, 5}, 23 ^ {0, 0, 6}, 24 ^ {2, 2, 3}, 25 ^ {1, 3, 3}, 

26 ^ {1, 2, 4}, 27 ^ {0, 3, 4}, 28 ^ {1, 1, 5}, 29 ^ {0, 2, 5}, 30 ^ {0, 1, 6}, 

31 ^ {0, 0, 7}, 32 ^ {2, 3, 3}, 33 ^ {2, 2, 4}, 34 ^ {1, 3, 4}, 35 ^ {0, 4, 4}, 
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36 ^{1,2, 5}, 37 ^{0,3, 5}, 38 ^{1,1, 6}, 39 ^{0,2, 6}, 40 ^{0,1, 7}, 

41 ^{0,0,8}, 42 ^{3, 3, 3}, 43 ^{2, 3, 4}, 44^ {1,4,4}, 45 ^{2, 2, 5}, 

46 ^{1,3, 5}, 47 ^{0,4, 5}, 48 ^{1,2, 6}, 49 ^{0,3, 6}, 50 ^{1,1, 7}, 

51 ^{0,2, 7}, 52 ^{0,1,8}, 53 ^{0,0, 9}. (A.l) 

Notice that we, in contrast to [27], start at index n = 1. 

B Shorthand notation for the bispectrum shapes 

In order to keep the different combinations of wavenumbers that contribute to the shape 
functions simple, we follow [28] and list the simplest terms that are consistent with the 
symmetries as follows: 



i 



K = Ki = fci + fes + fca, 



pq 



K — 



1 



A 



pqr 



(B.l) 
(B.2) 
(B.3) 



where A 
B.l 



— 1 + 5pq and A^^^ = Apq(Aqr + Spr). Spq is the Kronecker symbol. 
Conversion of K-terms 



Not ah of the ahowed K-terms in eq. (4.1) are independent of each other. Here we present 
the conversion for ah terms up to degree 6 
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expressed in terms of ad hoc chosen basic shapes (they happen to be highly correlated and 
thus in retrospective not an ideal choice). 
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(b) correlation 



Figure 29. (a) Plot of the sum of squared coefficients for the clipped local shape using different 
limits, (b) Plot of the correlation between the clipped shapes and the equilateral shape. 



C The clipping mechanism 

We checked the reliability of the clipping mechanism^^ proposed by Fergusson and Shellard 
[28] by applying different cuts, i.e. 

, ^\ < Xcut where Xcut ^ {0.030, 0.020, 0.010, 0.005, 0.015}, (C.l) 

to the local shape. The convergence of the expansion series in figure 29(a) indicates that the 
smaller the cut region the worse the convergence. However, the cut area should not be too 
large either in order to retain the physically relevant region in eq. (2.30), and to maintain 
a distinguishable local shape. The correlation with the equilateral shape, see figure 29(b), 
is a good indicator in that regard. Note that the partial correlation sums feature good 
convergence properties because the cut region does not contribute much to the correlation. 

The proposed cut Xcut = 0.015, which Fergusson and Shellard selected based on an 
investigation of the late-time behaviour, is satisfactory. 



D Late-time shapes 

We focused on the early-time bispectrum in this paper, given as 

N 



Bn{ki, A:2, k^) = , , ,^ S{kiM. ks)- (D.l) 



^^We chose the name "clipping mechanism" for the truncation of a shape to avoid misinterpretation with 
the truncation of an expansion series to some finite n. 
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However, for completeness, we would like to mention the relation to the CMB bispectrum 
(or late-time bispectrum) 

^ \mi 1712 msj 

which is a convolution of the early-time bispectrum with radiation transfer functions [12] 

X 9Th {h)9Ti2 {k2)9Th ih)jh {kix)ji^ {k2x)ji^ (ksx) , 

where gTi{k) is the temperature radiation transfer function for the /th multipole moment, 
jl{x) the corresponding Bessel function and we introduced the reduced bispectrum via the 
relation 

Bh,l,,ls=:Ql^[t,^,bi,i,i,. (D.4) 

Here, Ql^^J^^rn^ is a geometrical factor that depends only on the multipole moments and 
relates the isotropic reduced bispectrum to the full bispectrum. The convolution in eq. (D.3) 
leads to "acoustic peaks" in the late-time bispectrum, as discussed by Fergusson et al. in 
[29]. 

D.l /-domain expansion 

Fergusson et al. [27] identify the tetrahedal polynomials in the multipole domain, 

^1 ^ ^2 + ^3 for h ^ hi hi + cyclic perms., (D.5) 

h-^h^h^ 2n, n G N, 

where Zmax = 2000 is the maximal wavenumber observable by PLANCK. Motivated by the 
late-time estimator for /nl, an appropriate weight functions for this domain is 



which can be approximated by 
w{hihih) 



1 n.(2/^ + l)(2/ + |) /n^(2/-2/, + ^)e2^[i^g(2/)-iog(2m^^ 



2^2 ^^(2/_2/, + i) Y 21 + 1 + 1 ' (D.7) 

h + h + h = 2/, / E N. 

Note that there are some minor typographical errors in the version given by Shellard [27], 
eq. 46. Compare the differences in figure 30. 

Using the approximation, they observed that the 1-d polynomials qr{x) are almost 
the same as the A^-domain tetrahedral polynomials in figure 3 (see [27], figure 6, for a 
comparison). The next step is to use these polynomials to create the separable functions 
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500 1000 



Figure 30. Plot of the /-space weight function for the li = I2 = h = I case: 1), the original version, 
see eq. (D.6) (black, dashed), 2), our approximation in eq. (D.7) (green, thick), 3), the misprinted 
approximation in [27], eq. 46 (red), which, as we suppose, should be corrected to read as in eq. (D.7). 
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